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Abstract 

We perform classical molecular dynamics simulations to investigate the enhancement of the 
mass sensitivity and resonant frequency of graphene nanomechanical resonators that is achieved 
by driving them into the nonlinear oscillation regime. The mass sensitivity as measured by the 
resonant frequency shift is found to triple if the actuation energy is about 2.5 times the initial kinetic 
energy of the nanoresonator. The mechanism underlying the enhanced mass sensitivity is found to 
be the effective strain that is induced in the nanoresonator due to the nonlinear oscillations, where 
we obtain an analytic relationship between the induced effective strain and the actuation energy 
that is applied to the graphene nanoresonator. An important implication of this work is that there 
is no need for experimentalists to apply tensile strain to the resonators before actuation in order 
to enhance the mass sensitivity. Instead, enhanced mass sensitivity can be obtained by the far 
simpler technique of actuating nonlinear oscillations of an existing graphene nanoresonator. 
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I. INTRODUCTION 



Carbon nanotube and graphene-based nanomechanical resonators (GNMR) have at- 
tracted significant attention from the scientific community; see the recent reviews of Refs. [l- 
4|. Experiments have demonstrated that single-walled carbon nanotube nanomechanical 
resonators can serve as a mass sensor that is capable of detecting individual atoms or 
molecules^-, which is made possible by the high stiffness and low mass of carbon nan- 
otubes. Graphene also possesses high stiffness^— , but may prove to be a superior mass 
sensor than nanotubes due to its significantly larger surface area on which molecules can 
attach. 

The performance of GNMR mass sensors is closely related to its quality (Q)-factor, which 
reflects the energy that is lost during the mechanical oscillation of the resonator. In the first 
reported study of GNMRs, Bunch et al. observed very low (< 100) Q-factors for GNMRs 
working in the megahertz ranged. In a later experiment, the same group reported a dramatic 
increase of the Q-factor at lower temperatures, with the Q-factors reaching up to 9000 at 
10K for GNMRs produced from the chemical vapor deposition growth method^. Chen et al. 
also found that the Q-factor of GNMR increases with decreasing temperature, and reaches 
10 4 at 5 K— . More recently, Eichler et alM found that the Q-factors of GNMRs can reach 
values of 10 5 , which is very close to the theoretical Q-factor upper bound of about 10 6 , 
which has been predicted for GNMRs with all edges fixed coupled with idealized mechanical 
actuation^ 1 ^. 

Because a key sensing objective for GNMRs is to enable detection of individual molecules 
or atoms, it is critical to determine methods of enhancing their mass sensitivity. For exam- 
ple, Eichler et al. have shown that it is possible to increase the resonant frequency (and 
thus the mass sensitivity) of GNMRs by driving its mechanical oscillation into the nonlinear 
regime, which was explained using a continuum mechanics modeU 4 -. Using similar continuum 
models, several groups have theoretically studied the nonlinear effect on the mass sensitivity 
of different nanoresonators^— , with only a single, recent study related to graphene^. Be- 
sides the nonlinear resonance, mechanical strain is another possible method to improve the 
mass sensitivity of the GNMR, as the resonant frequency can be enhanced by mechanical 
tensiorjp£>2i~— . 

In this letter, we investigate the utility of inducing nonlinear oscillations to enhance 
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the mass sensitivity of GBMRs using classical molecular dynamics (MD) simulation. Our 
simulations show that the adsorption- induced frequency shift resulting from a single Au atom 
can be enhanced by increasing the actuation energy, when the actuation energy parameter 
a is below a critical value a c = 2.0875 ± 0.0125. For actuation energy above the critical 
value, the adsorbed Au atom exhibits significant diffusion, which reduces the Q-factor of 
the GNMR by more than three orders of magnitude and significantly degrades the mass 
sensitivity. Quite different from the widely used continuum mechanics explanation, we show 
that the nonlinear-induced frequency enhancement results from the 'effective strain' that is 
induced in the oscillating GNMR. We derive an analytic expression for the effective strain 
6 q = j 7r ' 2a m^L2 (here E^ is the kinetic energy, m the mass, uj the angular frequency, and L 
length of the resonator) that enables us to directly link the equivalence of applied mechanical 
tensile strain and the strain induced by nonlinear oscillations to the resonant frequencies of 
the GNMRs. 

II. SIMULATION DETAILS 

The graphene sample in our simulations has dimensions (L x ,L y ) = (34,38) A, and is 
composed of 504 carbon atoms. The atoms at the +x and — x ends of the GNMR are fixed, 
while periodic boundary conditions are applied in the y direction. The interactions of the 
carbon atoms are described by the Brenner (REBO-II) potential. For the cases where a 
single Au atom is adsorbed on the GNMR, the interaction between the Au atom and the 
GNMR is modeled by a Lennard- Jones potential with length parameter <r=2.9943 A and 
energy parameter e=0. 02936 eV— . The standard Newton equations of motion are integrated 
in time using the velocity Verlet algorithm with a time step of 1 fs. 

Our simulations are performed as follows. First, the Nose-Hoover— >21 thermostat is ap- 
plied to thermalize the system to a constant temperature of 4.2 K within the NVT ensemble, 
which is run for 10 5 MD steps. The mechanical oscillation of the resonator is then actuated 
by adding a velocity distribution to the system, which follows the morphology of the first 
flexural vibrational mode of graphene^. The imposed velocity distribution, or actuation en- 
ergy, is AE = aE®, where E% is the total kinetic energy in the GNMR after thermalization 
but just before its actuation and a is the actuation energy parameter. After the actuation 
energy is applied, the system is allowed to oscillate freely within the NVE ensemble for 2 19 
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MD steps. The data from the NVE ensemble is used to analyze the mechanical oscillation 
of the GNMR. 

III. RESULTS 

During the free vibration period, the energy in the GNMR switches between kinetic 
and potential energy. The frequency of the switching is 2f, with / being the frequency 
of mechanical oscillation. Fig. [TJ shows the resonant peaks that are obtained by taking a 
Fourier transformation of the time history of the kinetic energy; the resonant peaks are used 
to extract the resonant frequency / at liquid helium temperature (4.2 K) for GNMRs with 
different actuation energies a = 0.1, 0.5, 1.0, 1.5, and 2.0, where for reference, a = 1 means 
that the actuation energy equals the total kinetic energy in the GNMR after thermalization. 
It should be noted that the amplitude in Fig. [TJ is large, because of the long simulation time 
we have performed. We focus on the liquid helium temperature because this temperature is 
commonly utilized in experiments involving GNMRs (eg. 90 mK or 4.0 K in Ref. 

Fig. [2] shows the resonant frequency as a function of actuation energy a. Panel (a) 
shows that the resonant frequency is enhanced by increasing a. This enhancement is due 
to the increasingly nonlinear behavior of the GNMR induced by the increasing oscillation 
amplitude. Adsorption of a single Au atom on top of the GNMR causes a considerable 
reduction of the resonant frequency, due to the resulting increase in effective mass of the 
GNMR. This frequency shift is what is measured experimentally to detect the adsorbed mass, 
and it is seen in Panel (b) that the magnitude of the frequency shift increases by a factor 
of three for large actuation energies (a ~ 2.5) as compared to if the GNMR was actuated 
with a small a. Panel (b) shows that the frequency shift can be increased by applying larger 
actuation energies, which also results in increased mass sensitivity. However, as shown in 
the same figure, the frequency shift does not increase monotonically, and a decrease with 
increasing actuation energy is observed when a is large. 

The reduction in the frequency shift with increasing a is determined to be a result of 
diffusion of the adsorbed Au atom, which was previously observed to occur at GNMR 
temperatures exceeding about 30 K— . To confirm the diffusion mechanism, we calculate in 
Fig. El the mean free path for the Au atom as a function of a, where a sharp jump in the 
mean free path is observed at a critical value of actuation energy a c = 2.0875 ± 0.0125. For 
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a < ac c , the mean free path of the Au atom is around 0.5 A/ps with small fluctuations. In 
contrast, for a > a c , the mean free path increases to a value around 1.0 A/ps, which implies 
diffusion of the Au atom. 

The left top inset of Fig. |3]shows resonant curves for smaller actuation energy a = 0.1 and 
1.5. These smooth curves provide evidence that the mechanical oscillation of the GNMR is 
the only vibrational mode in the system. The right bottom inset of Fig. [3] shows a significant 
amount of diffusion-induced noise in the resonant curves of large actuation energy a = 4.0 
and 10.0. The noise corresponds to other vibrations that are induced by the diffusion of 
the Au atom. Fig. H] illustrates the trajectory history of the Au atom with a = 1.0, 2.075, 
2.1, and 4.0. The diffusion is clearly observed for actuation energy above the critical value 
a c . The thermal noise provides important energy damping channels for the mechanical 
oscillation of the GNMR. As a result, the Q-factor of the GNMR is greatly reduced as 
shown in Fig. |5j There is almost no energy dissipation in GNMR for actuation energy below 
a c , leading to extremely high Q-factors. For actuation energy above the critical value, the 
Q-factor is reduced by more than three orders of magnitude. The time history of the kinetic 
energy is shown in Fig. |5^b) for a = 2.075 and in Fig. [5](c) for a = 2.1, from which the 
Q-factor has been extracted^. 

Fig. [6] shows the Q-factor vs. the actuation energy parameter a for graphene nanome- 
chanical resonators with three different adsorbates. In these three systems, the adsorbate 
mass m and the Lennard- Jones interaction parameter e are (m, e), (m/2, e), and (m, e/2). 
For the first system (m, e), the mass of the adsorbate is m = 197 and the Lennard- Jones 
potential parameter is e = 0.02936 eV, which are the actual values for Au adsorbate. For 
the second system (m/2, e), the mass of the adsorbate is reduced by half to be m = 99 
while e (the interaction strength between the adsorbate and graphene) remains unchanged. 
In the third system (m, e/2), the mass of the adsorbate is unchanged at m = 197 while the 
interaction strength e between the adsorbate and graphene is reduced by half. 

In all three systems, there is a sharp decrease in the Q-factor at different values of actu- 
ation energy parameter a, which results from the diffusion of the adsorbate. These results 
also follow physical intuition. For example, diffusion happens at smaller a for adsorbates 
that have a weaker bonding strength with graphene. In contrast, diffusion happens at larger 
a for adsorbates with less mass, since those atoms have a smaller kinetic energy at the same 
temperature as larger mass adsorbates, and therefore more energy via the nonlinear actu- 
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ation parameter a is needed to induce diffusion for the smaller mass atoms. These results 
may serve as useful guidelines for experimentalists to verify our theoretical predictions, as 
different types of adsorbates (of different mass or/and different interaction strengths with 
graphene) are usually observed experimentally 

IV. DISCUSSION 

The above discussion has established that the mass sensitivity of the GNMR can be 
enhanced by driving the mechanical oscillations into the nonlinear regime using a larger 
actuation energy. The remainder of this article will be devoted to explaining the mechanism 
that enables the nonlinear-induced resonant frequency, and thus mass sensitivity enhance- 
ment. 

As illustrated in Fig. the GNMR is initially flat (gray points on the horizontal line) with 
length L. After exciting the mechanical oscillations via the actuation energy AE = aE®, 
the GNMR oscillates with amplitude A = v/f^, where to = 2?rf is the angular frequency 
of the mechanical oscillation and m is the total mass of the system. In the derivation of the 
oscillation amplitude, the thermal vibrations at 4.2 K have been ignored due to the very 
low temperature conditions. We focus on a particular point in the GNMR, i.e. the atom at 
the midpoint, which oscillates as u = A sin ut. The mean oscillation amplitude of this point 
is V < u 2 > = A/y/2, and so this point can effectively be regarded as a stationary point 
located at A/y/2. Similarly, the shape of the oscillating GNMR is equivalent to a stationary 
sine function, but with amplitude A/y2. This effective shape is shown as big red points in 
Fig. [7J where the length of the effective shape is 

S= ^Ly/^TS (K(k)), (1) 
7T 

where a = ttA/\/2L, k = a 2 /(a 2 + 1) and K{k) is the complete elliptic integral of the first 
kind. 

In analyzing the sine-shaped form of the oscillating GNMR (red dots) in Fig. [7J our 
key insight is that the oscillating GNMR is under an 'effective strain' as compared to the 
original, flat GNMR (gray points on the horizontal line). This effective strain is denoted by 
e a , where the subscript indicates that the effective strain is induced by the actuation energy 
a, and where e a = (S — L)/L. Approximating the elliptic integral up to second order, we 
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can derive an analytic expression for the oscillation-induced effective strain: 

3 E° 

e » = -A^ a TF2- W 

4 mu 2 L z 

Substituting all structural and physical parameters into the expression, we get a concise rela- 
tionship between the effective strain and the actuation energy that is applied to the GNMR: 
e a ~ 0.216% x a. The effective strain for the Au adsorbed GNMR is obtained analogously 
As previously discussed, the resonant frequency of the GNMR can be enhanced by applying 
mechanical tension due to the enhanced stiffening of the structure that results^^— . In this 
sense, we propose that the enhancement of the resonant frequency that results from driving 
the GNMR into the nonlinear oscillation regime is actually due to the effective strain that 
is induced in the oscillating GNMR. 

In Fig. El we show that the effective strain is indeed responsible for the enhancement of 
the resonant frequency with increasing a. We present the resonant frequency for the GNMR 
in (a) and the GNMR with adsorption of a single Au atom in (b). The circles represent 
the resonant frequency of the GNMR due to different actuation energies a, i.e f(a), while 
the dashed lines (blue online) are the resonant frequencies f(e a ) for the GNMR that was 
pre-stretched with a tensile strain e a before actuation, where the strain e a for actuation 
energy a is calculated by Eq. (J2J). A very small actuation energy of a = 0.001 is used to 
actuate the stretched GNMR such only the strain effect is important in the GNMR under 
mechanical tension. 

Both Figs. IU(a) and (b) show good agreement between the resonant frequencies obtained 
from increasing the actuation energy a and by applying tensile strain e a . This agreement 
verifies that the effective strain is the cause for the increase of the resonant frequency that we 
have previously observed by increasing a. The agreement is particularly good for a < 1.0, 
though there is an increasing discrepancy for larger a. The solid lines (red online) show that 
this discrepancy is greatly improved by considering a high-order correction in the square of 
the oscillation amplitude A 2 . This high-order correction is demonstrated in Fig. [8](c), where 
u z on the y-axis of Fig. [8](c) is the oscillation of the midpoint of the GNMR obtained from 
the MD simulation, and where Fig. [8fc) is plotted using a log-log scale. < u 2 z > deviates 
from linear behavior because of the phonon-phonon scattering. Although the temperature is 
quite low (4.2 K), the phonon-phonon scattering is still important for the first bending mode 
(i.e the mode of mechanical oscillation) because this mode is in a highly non-equilibrium 



state. In other words, because the first bending mode has been mechanically actuated with 
very large amplitude, it is driven far away from its thermal equilibrium state at 4.2 K. 
Our MD simulations gives < u 2 >oc a 0,79 for the pure GNMR. However, we have assumed 
< u 2 >oc a 10 in our derivation of the effective strain. The effective strain e a oc< u 2 >, which 
is why we have e a oc a. Therefore, the effective strain must be modified by: e' a = e a /a ' 21 , 
where e Q is the effective strain without correction. A similar correction factor of a ' 28 can 
be done for the effective strain of the GNMR with a single adsorbed Au atom. 

We note that Fig. [S^c) shows a small oscillation amplitude ( 0.8 A) for the system sim- 
ulated in this work, due to its small size. The amplitude U actually depends on the length 
L of the graphene as U = 2.43 x A, where the frequency u is length-dependent. For 
the system we have simulated, this formula gives U ~ 0.8 A at a = 1.0, which is in good 
agreement with MD simulation results in Fig. [S^c). Considering the flexural property of 
the bending mode in graphene, i.e u oc 1/L 2 , we can obtain the oscillation amplitude of 
an arbitrary system length as U = Uq(L / Lq) 2 , where Lq = 34 A and Uq = 0.8 A for our 
simulated system at a = 1.0. For instance, the mechanical oscillation amplitude becomes 
67 nm if the system length is 100 nm. If the system length is 1 /xm, then the mechanical os- 
cillation amplitude is very large, i.e. 6.7 /xm. These results show that a = 1 is large enough 
to induce oscillations that can be clearly distinguished from random thermal fluctuations in 
experiments. 

Previously, we have established that the effective strain mechanism can explain the en- 
hancement of the resonant frequency that occurs by increasing a. However, we also needed 
to perform a high-order correction for large actuation energies to account for the phonon- 
phonon scattering that is still present even at the low (4.2 K) temperature due to the highly 
non-equilibrium behavior resulting from the large actuation energy. However, if the temper- 
ature is reduced to K, there is no phonon-phonon scattering for the first bending mode 
because this mode cannot decay without the assistance of other vibrations according to 
symmetry selection rules^. As a result, at K, the enhancement of the resonant frequency 
by applying tensile strain before actuation should agree with that obtained by increasing 
the actuation energy, without any high-order correction. 

Indeed, we observe this result in Fig. [8](d). We emphasize that for these two simulations 
of the GNMR, the GNMR exists initially in an energy minimized configuration at K, and 
the applied actuation energy is AE = aE®. Because the kinetic energy E® is at K, we 
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use E% = Ek(4:.2K), where Ek(4:.2K) is the total kinetic energy at 4.2 K, for both of the two 
simulations shown in Fig. [8] (d). The solid line is the kinetic energy for the GNMR that is 
actuated by the large actuation energy parameter a = 2.0. The dashed line (blue online) is 
the kinetic energy of the GNMR that is pre-stretched in tension by 0.4% before actuation, 
where 0.4% is chosen as it is also the effective strain that results from the actuation energy 
a = 2.0. The pre-stretched GNMR is actuated with a very small actuation energy a = 0.001 
so that the effect from a can be ignored and only strain is important in this case. We observe 
exactly the same resonant frequency in the mechanical oscillation in these two very different 
situations. From the Fourier transformation of the time history of thees two kinetic energy, 
we get the same resonant frequency of 270 GHz, which further establishes the equivalent 
effect of applied tensile strain and different actuation energies on the resonant frequencies 
of GNMRs. 

V. CONCLUSION 

In conclusion, we used classical MD simulations to demonstrate that a significant en- 
hancement of the mass sensitivity of graphene nanomechanical resonators can be achieved 
by driving them into the nonlinear oscillation regime. The enabling mechanism was deter- 
mined to be the effective strain that is induced in the nanoresonators due to the nonlinear 
oscillations. A simple analytic expression relating the effective strain to the actuation energy 
was obtained, and shown to be quite accurate for moderate actuation energies. The key im- 
plication is that it should be possible for experimentalists to directly incorporate the present 
findings to enhance the sensitivity of graphene-based mass sensors simply by actuating the 
graphene nanomechanical resonators into the nonlinear oscillation regime. 
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FIG. 1: (Color online) Resonant frequencies for GNMRs. From left to right are pure GNMR at 
4.2 K with different actuation energy parameters a = 0.1, 0.5, 1.0, 1.5, and 2.0. 
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FIG. 2: (Color online) The resonant frequency versus actuation energy parameter a. (a) The reso- 
nant frequency for pure GNMR (circle points) and GNMR adsorbed by one Au atom (pentagonal 
points) . (b) . The magnitude of frequency shift induced by the adsorption of a single Au atom for 
different actuation energy parameters a. 
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FIG. 3: (Color online) The mean free path of the Au atom as a function of the actuation energy 
parameter a. A sharp jump in the mean free path occurs at a c = 2.0875 ± 0.0125, implying the 
diffusion of the Au atom for a > a c . Left top inset shows smooth resonant curves for smaller 
actuation energy parameters a = 0.1 and 1.5. Right bottom inset shows significant diffusion- 
induced noise in the resonant curves for large actuation energy parameters a = 4.0 and 10.0. 
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FIG. 4: The trajectory history of the adsorbed Au atom for actuation energy a = 1.0, 2.075, 2.1, 
and 4.0 from left to right. The first two panels show a localized vibration of the Au atom, while 
the last two panels clearly show diffusion of the Au atom. 
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FIG. 5: (Color online) The effect of the diffusion of the Au atom on the Q-factor of the GNMR. 
(a). The Q-factors of Au adsorbed GNMR remains extremely high value for a < a c , but becomes 
three orders of magnitude smaller for a greater than the critical value a c . The time history of the 
total kinetic energy is shown in (b) for a = 2.075 and in (c) for a = 2.1, which are very close to 
a c . The energy decay in (c) is caused by the diffusion of the adsorbed Au atom shown in Fig. HI 
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FIG. 6: (Color online) Quality factors versus the actuation energy parameter a for graphene 
nanomechanical resonators adsorbed by a single atom with mass and Lennard- Jones potential 
parameter as (m, e), (m/2, e), and (m, e/2). 




FIG. 7: (Color online) Geometry of the GNMR. Before actuation, the GNMR lies in the horizontal 
plane as shown by the small gray points on the horizontal line. After mechanical actuation, the 
GNMR oscillates with amplitude A, which is determined by the actuation energy (see text for more 
details). The oscillating GNMR is regarded as an effective stationary shape with amplitude A/\/2 
(big red points). The difference between the length of the effective shape (S) and the length of the 
initial structure (L) yields an effective strain e a = (S — L)/L. 
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FIG. 8: (Color online) Resonant frequencies for pure GNMR in (a) and GNMR with adsorption of 
an Au atom in (b). The resonant frequency can be enhanced by increasing actuation energy a, i.e 
f(a) (circle points); or by applying effective strain, i.e f(e a ) (dashed blue lines); or by applying a 
modified effective strain, i.e f(e' a ) (solid red lines), (c) The square of the mean vibration amplitude 
for the atom in the center of the GNMR from the MD simulations; Lnx and Lny signify log e x and 
log e y, respectively, (d) The time history of the kinetic energy for GNMR with actuation energy 
parameter a = 2.0 or for the GNMR under mechanical tension e a = 0.4%, both at K. 
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